\(\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} (B \sec (c+d x)+C \sec ^2(c+d x)) \, dx\) [376]

Optimal result

Integrand size = 40, antiderivative size = 175 \[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\frac {64 a^3 (15 B+13 C) \tan (c+d x)}{315 d \sqrt {a+a \sec (c+d x)}}+\frac {16 a^2 (15 B+13 C) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a (15 B+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}+\frac {2 (9 B-2 C) (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d} \]

[Out]

2/105*a*(15*B+13*C)*(a+a*sec(d*x+c))^(3/2)*tan(d*x+c)/d+2/63*(9*B-2*C)*(a+a*sec(d*x+c))^(5/2)*tan(d*x+c)/d+2/9
*C*(a+a*sec(d*x+c))^(7/2)*tan(d*x+c)/a/d+64/315*a^3*(15*B+13*C)*tan(d*x+c)/d/(a+a*sec(d*x+c))^(1/2)+16/315*a^2
*(15*B+13*C)*(a+a*sec(d*x+c))^(1/2)*tan(d*x+c)/d

Rubi [A] (verified)

Time = 0.50 (sec) , antiderivative size = 175, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.125, Rules used = {4157, 4095, 4086, 3878, 3877} \[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\frac {64 a^3 (15 B+13 C) \tan (c+d x)}{315 d \sqrt {a \sec (c+d x)+a}}+\frac {16 a^2 (15 B+13 C) \tan (c+d x) \sqrt {a \sec (c+d x)+a}}{315 d}+\frac {2 (9 B-2 C) \tan (c+d x) (a \sec (c+d x)+a)^{5/2}}{63 d}+\frac {2 a (15 B+13 C) \tan (c+d x) (a \sec (c+d x)+a)^{3/2}}{105 d}+\frac {2 C \tan (c+d x) (a \sec (c+d x)+a)^{7/2}}{9 a d} \]

[In]

Int[Sec[c + d*x]*(a + a*Sec[c + d*x])^(5/2)*(B*Sec[c + d*x] + C*Sec[c + d*x]^2),x]

[Out]

(64*a^3*(15*B + 13*C)*Tan[c + d*x])/(315*d*Sqrt[a + a*Sec[c + d*x]]) + (16*a^2*(15*B + 13*C)*Sqrt[a + a*Sec[c
+ d*x]]*Tan[c + d*x])/(315*d) + (2*a*(15*B + 13*C)*(a + a*Sec[c + d*x])^(3/2)*Tan[c + d*x])/(105*d) + (2*(9*B
- 2*C)*(a + a*Sec[c + d*x])^(5/2)*Tan[c + d*x])/(63*d) + (2*C*(a + a*Sec[c + d*x])^(7/2)*Tan[c + d*x])/(9*a*d)

Rule 3877

Int[csc[(e_.) + (f_.)*(x_)]*Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Simp[-2*b*(Cot[e + f*x]/(
f*Sqrt[a + b*Csc[e + f*x]])), x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0]

Rule 3878

Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> Simp[(-b)*Cot[e + f*x]*(
(a + b*Csc[e + f*x])^(m - 1)/(f*m)), x] + Dist[a*((2*m - 1)/m), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^(m - 1),
x], x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0] && GtQ[m, 1/2] && IntegerQ[2*m]

Rule 4086

Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*(B_.) + (A_))
, x_Symbol] :> Simp[(-B)*Cot[e + f*x]*((a + b*Csc[e + f*x])^m/(f*(m + 1))), x] + Dist[(a*B*m + A*b*(m + 1))/(b
*(m + 1)), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m, x], x] /; FreeQ[{a, b, A, B, e, f, m}, x] && NeQ[A*b - a*B
, 0] && EqQ[a^2 - b^2, 0] && NeQ[a*B*m + A*b*(m + 1), 0] &&  !LtQ[m, -2^(-1)]

Rule 4095

Int[csc[(e_.) + (f_.)*(x_)]^2*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*(B_.) + (A_
)), x_Symbol] :> Simp[(-B)*Cot[e + f*x]*((a + b*Csc[e + f*x])^(m + 1)/(b*f*(m + 2))), x] + Dist[1/(b*(m + 2)),
Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m*Simp[b*B*(m + 1) + (A*b*(m + 2) - a*B)*Csc[e + f*x], x], x], x] /; Fr
eeQ[{a, b, e, f, A, B, m}, x] && NeQ[A*b - a*B, 0] &&  !LtQ[m, -1]

Rule 4157

Int[((a_.) + csc[(e_.) + (f_.)*(x_)]*(b_.))^(m_.)*((A_.) + csc[(e_.) + (f_.)*(x_)]*(B_.) + csc[(e_.) + (f_.)*(
x_)]^2*(C_.))*((c_.) + csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_.), x_Symbol] :> Dist[1/b^2, Int[(a + b*Csc[e + f*x])
^(m + 1)*(c + d*Csc[e + f*x])^n*(b*B - a*C + b*C*Csc[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C, m,
n}, x] && EqQ[A*b^2 - a*b*B + a^2*C, 0]

Rubi steps \begin{align*} \text {integral}& = \int \sec ^2(c+d x) (a+a \sec (c+d x))^{5/2} (B+C \sec (c+d x)) \, dx \\ & = \frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {2 \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (\frac {7 a C}{2}+\frac {1}{2} a (9 B-2 C) \sec (c+d x)\right ) \, dx}{9 a} \\ & = \frac {2 (9 B-2 C) (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {1}{21} (15 B+13 C) \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \, dx \\ & = \frac {2 a (15 B+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}+\frac {2 (9 B-2 C) (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {1}{105} (8 a (15 B+13 C)) \int \sec (c+d x) (a+a \sec (c+d x))^{3/2} \, dx \\ & = \frac {16 a^2 (15 B+13 C) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a (15 B+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}+\frac {2 (9 B-2 C) (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {1}{315} \left (32 a^2 (15 B+13 C)\right ) \int \sec (c+d x) \sqrt {a+a \sec (c+d x)} \, dx \\ & = \frac {64 a^3 (15 B+13 C) \tan (c+d x)}{315 d \sqrt {a+a \sec (c+d x)}}+\frac {16 a^2 (15 B+13 C) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a (15 B+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}+\frac {2 (9 B-2 C) (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.81 (sec) , antiderivative size = 96, normalized size of antiderivative = 0.55 \[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\frac {2 a^3 \left (690 B+584 C+(345 B+292 C) \sec (c+d x)+3 (60 B+73 C) \sec ^2(c+d x)+5 (9 B+26 C) \sec ^3(c+d x)+35 C \sec ^4(c+d x)\right ) \tan (c+d x)}{315 d \sqrt {a (1+\sec (c+d x))}} \]

[In]

Integrate[Sec[c + d*x]*(a + a*Sec[c + d*x])^(5/2)*(B*Sec[c + d*x] + C*Sec[c + d*x]^2),x]

[Out]

(2*a^3*(690*B + 584*C + (345*B + 292*C)*Sec[c + d*x] + 3*(60*B + 73*C)*Sec[c + d*x]^2 + 5*(9*B + 26*C)*Sec[c +
d*x]^3 + 35*C*Sec[c + d*x]^4)*Tan[c + d*x])/(315*d*Sqrt[a*(1 + Sec[c + d*x])])

Maple [A] (verified)

Time = 72.18 (sec) , antiderivative size = 133, normalized size of antiderivative = 0.76

method	result	size
default	\(\frac {2 a^{2} \left (690 B \cos \left (d x +c \right )^{4}+584 C \cos \left (d x +c \right )^{4}+345 B \cos \left (d x +c \right )^{3}+292 C \cos \left (d x +c \right )^{3}+180 B \cos \left (d x +c \right )^{2}+219 C \cos \left (d x +c \right )^{2}+45 B \cos \left (d x +c \right )+130 C \cos \left (d x +c \right )+35 C \right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{3}}{315 d \left (\cos \left (d x +c \right )+1\right )}\)	\(133\)
parts	\(\frac {2 B \,a^{2} \left (46 \cos \left (d x +c \right )^{3}+23 \cos \left (d x +c \right )^{2}+12 \cos \left (d x +c \right )+3\right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{2}}{21 d \left (\cos \left (d x +c \right )+1\right )}+\frac {2 C \,a^{2} \left (584 \cos \left (d x +c \right )^{4}+292 \cos \left (d x +c \right )^{3}+219 \cos \left (d x +c \right )^{2}+130 \cos \left (d x +c \right )+35\right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{3}}{315 d \left (\cos \left (d x +c \right )+1\right )}\)	\(162\)

[In]

int(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(B*sec(d*x+c)+C*sec(d*x+c)^2),x,method=_RETURNVERBOSE)

[Out]

2/315*a^2/d*(690*B*cos(d*x+c)^4+584*C*cos(d*x+c)^4+345*B*cos(d*x+c)^3+292*C*cos(d*x+c)^3+180*B*cos(d*x+c)^2+21
9*C*cos(d*x+c)^2+45*B*cos(d*x+c)+130*C*cos(d*x+c)+35*C)*(a*(1+sec(d*x+c)))^(1/2)/(cos(d*x+c)+1)*tan(d*x+c)*sec
(d*x+c)^3

Fricas [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 136, normalized size of antiderivative = 0.78 \[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\frac {2 \, {\left (2 \, {\left (345 \, B + 292 \, C\right )} a^{2} \cos \left (d x + c\right )^{4} + {\left (345 \, B + 292 \, C\right )} a^{2} \cos \left (d x + c\right )^{3} + 3 \, {\left (60 \, B + 73 \, C\right )} a^{2} \cos \left (d x + c\right )^{2} + 5 \, {\left (9 \, B + 26 \, C\right )} a^{2} \cos \left (d x + c\right ) + 35 \, C a^{2}\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sin \left (d x + c\right )}{315 \, {\left (d \cos \left (d x + c\right )^{5} + d \cos \left (d x + c\right )^{4}\right )}} \]

[In]

integrate(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(B*sec(d*x+c)+C*sec(d*x+c)^2),x, algorithm="fricas")

[Out]

2/315*(2*(345*B + 292*C)*a^2*cos(d*x + c)^4 + (345*B + 292*C)*a^2*cos(d*x + c)^3 + 3*(60*B + 73*C)*a^2*cos(d*x
+ c)^2 + 5*(9*B + 26*C)*a^2*cos(d*x + c) + 35*C*a^2)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sin(d*x + c)/(d*
cos(d*x + c)^5 + d*cos(d*x + c)^4)

Sympy [F]

\[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\int \left (a \left (\sec {\left (c + d x \right )} + 1\right )\right )^{\frac {5}{2}} \left (B + C \sec {\left (c + d x \right )}\right ) \sec ^{2}{\left (c + d x \right )}\, dx \]

[In]

integrate(sec(d*x+c)*(a+a*sec(d*x+c))**(5/2)*(B*sec(d*x+c)+C*sec(d*x+c)**2),x)

[Out]

Integral((a*(sec(c + d*x) + 1))**(5/2)*(B + C*sec(c + d*x))*sec(c + d*x)**2, x)

Maxima [F]

\[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\int { {\left (C \sec \left (d x + c\right )^{2} + B \sec \left (d x + c\right )\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {5}{2}} \sec \left (d x + c\right ) \,d x } \]

[In]

integrate(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(B*sec(d*x+c)+C*sec(d*x+c)^2),x, algorithm="maxima")

[Out]

4/315*(315*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)*((B*a^2*d*cos(2*d*x + 2*c)
^4 + B*a^2*d*sin(2*d*x + 2*c)^4 + 4*B*a^2*d*cos(2*d*x + 2*c)^3 + 6*B*a^2*d*cos(2*d*x + 2*c)^2 + 4*B*a^2*d*cos(
2*d*x + 2*c) + B*a^2*d + 2*(B*a^2*d*cos(2*d*x + 2*c)^2 + 2*B*a^2*d*cos(2*d*x + 2*c) + B*a^2*d)*sin(2*d*x + 2*c
)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)
*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*
c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2
*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*
x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*
d*x + 2*c))*sin(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x
+ 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*
sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c
)*sin(2*d*x + 2*c))*cos(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c)
+ 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x +
8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x
+ 2*c)^2)*sin(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x +
2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6
*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*co
s(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(
cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 +
sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2
+ sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos
(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2
*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(
2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*
sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (co
s(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos
(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(
2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2
+ (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2
+ 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c
)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2
*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*
x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)
^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2
*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^
2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6
*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c
)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)
^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x +
2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)), x) + ((11*B + 10*C)*a^2*d*cos(2*d*x + 2*c)^4 +
(11*B + 10*C)*a^2*d*sin(2*d*x + 2*c)^4 + 4*(11*B + 10*C)*a^2*d*cos(2*d*x + 2*c)^3 + 6*(11*B + 10*C)*a^2*d*cos(
2*d*x + 2*c)^2 + 4*(11*B + 10*C)*a^2*d*cos(2*d*x + 2*c) + (11*B + 10*C)*a^2*d + 2*((11*B + 10*C)*a^2*d*cos(2*d
*x + 2*c)^2 + 2*(11*B + 10*C)*a^2*d*cos(2*d*x + 2*c) + (11*B + 10*C)*a^2*d)*sin(2*d*x + 2*c)^2)*integrate((((c
os(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) +
cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x +
4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x
+ 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*
x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(7/2
*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((co
s(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) -
cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))
*cos(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c
)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2
*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(7/2*
arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/(((cos
(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(
8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(
cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 +
(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2
+ sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*
c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^
2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + c
os(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(
4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x +
2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*
(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) +
3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*s
in(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(
2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)
^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^
2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2
*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x +
2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*
x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos
(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d
*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*
c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)
+ 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d
*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x +
2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x +
2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*
c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*
c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))
*sin(4*d*x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x
+ 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)), x) + 5*((3*B + 4*C)*a^2*d*cos(2*d*x + 2*c)^4 + (3*B + 4*C)*a^2*d*s
in(2*d*x + 2*c)^4 + 4*(3*B + 4*C)*a^2*d*cos(2*d*x + 2*c)^3 + 6*(3*B + 4*C)*a^2*d*cos(2*d*x + 2*c)^2 + 4*(3*B +
4*C)*a^2*d*cos(2*d*x + 2*c) + (3*B + 4*C)*a^2*d + 2*((3*B + 4*C)*a^2*d*cos(2*d*x + 2*c)^2 + 2*(3*B + 4*C)*a^2
*d*cos(2*d*x + 2*c) + (3*B + 4*C)*a^2*d)*sin(2*d*x + 2*c)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x + 2*c) +
3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8
*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x +
2*c)^2)*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d
*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6
*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*
x + 2*c))))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3
*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) -
3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(5/2*arctan2(sin(2*d*x + 2*c),
cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x +
4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x
+ 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x
+ 2*c))))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^
4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c
)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*
c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2
*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*
x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(
2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d
*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*
x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2
+ 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2
+ cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*
d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x +
2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2
*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) +
6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos
(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*
x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d
*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(
2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*s
in(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 +
(2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*
sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3
*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 +
cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*
x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x
+ 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2
+ cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + si
n(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*c
os(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*
d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(
4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x
+ 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2*arctan2(
sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)
^(1/4)), x) + ((5*B + 2*C)*a^2*d*cos(2*d*x + 2*c)^4 + (5*B + 2*C)*a^2*d*sin(2*d*x + 2*c)^4 + 4*(5*B + 2*C)*a^2
*d*cos(2*d*x + 2*c)^3 + 6*(5*B + 2*C)*a^2*d*cos(2*d*x + 2*c)^2 + 4*(5*B + 2*C)*a^2*d*cos(2*d*x + 2*c) + (5*B +
2*C)*a^2*d + 2*((5*B + 2*C)*a^2*d*cos(2*d*x + 2*c)^2 + 2*(5*B + 2*C)*a^2*d*cos(2*d*x + 2*c) + (5*B + 2*C)*a^2
*d)*sin(2*d*x + 2*c)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) +
3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x +
6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(3/2*arctan2(sin(2*d*x +
2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*
d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(
4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin(2*d*x
+ 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) +
3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c)
- 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c
)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2
*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x
+ 2*c) + sin(2*d*x + 2*c)^2)*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(2*d*x
+ 2*c), cos(2*d*x + 2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*
x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*
d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8
*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(c
os(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(
cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*
d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c)
+ 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(
4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin
(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*
x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*
x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*
d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*c
os(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)
*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2
+ 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2
+ 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*
c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x +
2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2
*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*
c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) +
6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*c
os(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d
*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c)
+ cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)
+ (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 +
3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*c
os(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) +
1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)), x))*sqrt(a) - (3*(35*(5*B +
2*C)*a^2*sin(6*d*x + 6*c) + 21*(25*B + 24*C)*a^2*sin(4*d*x + 4*c) + 3*(155*B + 146*C)*a^2*sin(2*d*x + 2*c))*co
s(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - (105*(5*B + 2*C)*a^2*cos(6*d*x + 6*c) + 63*(25*B + 24
*C)*a^2*cos(4*d*x + 4*c) + 9*(155*B + 146*C)*a^2*cos(2*d*x + 2*c) + (345*B + 292*C)*a^2)*sin(9/2*arctan2(sin(2
*d*x + 2*c), cos(2*d*x + 2*c) + 1)))*sqrt(a))/((d*cos(2*d*x + 2*c)^4 + d*sin(2*d*x + 2*c)^4 + 4*d*cos(2*d*x +
2*c)^3 + 6*d*cos(2*d*x + 2*c)^2 + 2*(d*cos(2*d*x + 2*c)^2 + 2*d*cos(2*d*x + 2*c) + d)*sin(2*d*x + 2*c)^2 + 4*d
*cos(2*d*x + 2*c) + d)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4))

Giac [F]

\[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\int { {\left (C \sec \left (d x + c\right )^{2} + B \sec \left (d x + c\right )\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {5}{2}} \sec \left (d x + c\right ) \,d x } \]

[In]

integrate(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(B*sec(d*x+c)+C*sec(d*x+c)^2),x, algorithm="giac")

[Out]

sage0*x

Mupad [B] (verification not implemented)

Time = 28.03 (sec) , antiderivative size = 723, normalized size of antiderivative = 4.13 \[ \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx=\frac {\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (\frac {B\,a^2\,4{}\mathrm {i}}{3\,d}-\frac {a^2\,\left (60\,B+73\,C\right )\,8{}\mathrm {i}}{315\,d}\right )+\frac {a^2\,\left (5\,B+2\,C\right )\,4{}\mathrm {i}}{3\,d}\right )\,\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}+\frac {\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (-\frac {B\,a^2\,4{}\mathrm {i}}{5\,d}+\frac {a^2\,\left (3\,B+4\,C\right )\,16{}\mathrm {i}}{105\,d}+\frac {a^2\,\left (9\,B+10\,C\right )\,4{}\mathrm {i}}{5\,d}\right )-\frac {a^2\,\left (5\,B+2\,C\right )\,4{}\mathrm {i}}{5\,d}+\frac {a^2\,\left (5\,B+16\,C\right )\,4{}\mathrm {i}}{5\,d}\right )}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,{\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}^2}+\frac {\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (\frac {B\,a^2\,4{}\mathrm {i}}{7\,d}+\frac {C\,a^2\,32{}\mathrm {i}}{63\,d}+\frac {a^2\,\left (B+2\,C\right )\,20{}\mathrm {i}}{7\,d}-\frac {a^2\,\left (B+C\right )\,40{}\mathrm {i}}{7\,d}\right )+\frac {a^2\,\left (B-8\,C\right )\,4{}\mathrm {i}}{7\,d}+\frac {a^2\,\left (5\,B+2\,C\right )\,4{}\mathrm {i}}{7\,d}-\frac {a^2\,\left (5\,B+9\,C\right )\,8{}\mathrm {i}}{7\,d}\right )}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,{\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}^3}-\frac {\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (\frac {B\,a^2\,4{}\mathrm {i}}{9\,d}+\frac {a^2\,\left (3\,B+4\,C\right )\,20{}\mathrm {i}}{9\,d}-\frac {a^2\,\left (5\,B+2\,C\right )\,4{}\mathrm {i}}{9\,d}-\frac {a^2\,\left (11\,B+10\,C\right )\,4{}\mathrm {i}}{9\,d}\right )-\frac {B\,a^2\,4{}\mathrm {i}}{9\,d}-\frac {a^2\,\left (3\,B+4\,C\right )\,20{}\mathrm {i}}{9\,d}+\frac {a^2\,\left (5\,B+2\,C\right )\,4{}\mathrm {i}}{9\,d}+\frac {a^2\,\left (11\,B+10\,C\right )\,4{}\mathrm {i}}{9\,d}\right )}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,{\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}^4}-\frac {a^2\,{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left (345\,B+292\,C\right )\,4{}\mathrm {i}}{315\,d\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )} \]

[In]

int(((B/cos(c + d*x) + C/cos(c + d*x)^2)*(a + a/cos(c + d*x))^(5/2))/cos(c + d*x),x)

[Out]

((exp(c*1i + d*x*1i)*((B*a^2*4i)/(3*d) - (a^2*(60*B + 73*C)*8i)/(315*d)) + (a^2*(5*B + 2*C)*4i)/(3*d))*(a + a/
(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)) +
((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(exp(c*1i + d*x*1i)*((a^2*(3*B + 4*C)*16i)/(105
*d) - (B*a^2*4i)/(5*d) + (a^2*(9*B + 10*C)*4i)/(5*d)) - (a^2*(5*B + 2*C)*4i)/(5*d) + (a^2*(5*B + 16*C)*4i)/(5*
d)))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)^2) + ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*
1i)/2))^(1/2)*(exp(c*1i + d*x*1i)*((B*a^2*4i)/(7*d) + (C*a^2*32i)/(63*d) + (a^2*(B + 2*C)*20i)/(7*d) - (a^2*(B
+ C)*40i)/(7*d)) + (a^2*(B - 8*C)*4i)/(7*d) + (a^2*(5*B + 2*C)*4i)/(7*d) - (a^2*(5*B + 9*C)*8i)/(7*d)))/((exp
(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)^3) - ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1
/2)*(exp(c*1i + d*x*1i)*((B*a^2*4i)/(9*d) + (a^2*(3*B + 4*C)*20i)/(9*d) - (a^2*(5*B + 2*C)*4i)/(9*d) - (a^2*(1
1*B + 10*C)*4i)/(9*d)) - (B*a^2*4i)/(9*d) - (a^2*(3*B + 4*C)*20i)/(9*d) + (a^2*(5*B + 2*C)*4i)/(9*d) + (a^2*(1
1*B + 10*C)*4i)/(9*d)))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)^4) - (a^2*exp(c*1i + d*x*1i)*(a + a
/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(345*B + 292*C)*4i)/(315*d*(exp(c*1i + d*x*1i) + 1))